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Tj" ■ PERTURBED RISK PROCESSES 1 

o 

^ ■ By Miljenko Huzak, 2 Mihael Perman, Hrvoje Sikic 3 

and zoran vondracek 4 

University of Zagreb, University of Ljubljana and University of Zagreb 

We study a general perturbed risk process with cumulative claims 
modelled by a subordinator with finite expectation, with the pertur- 
|"y/| bation being a spectrally negative Levy process with zero expecta- 

tion. We derive a Pollaczek-Hinchin type formula for the survival 
probability of that risk process, and give an interpretation of the for- 
| mula based on the decomposition of the dual risk process at modified 

, ladder epochs. 

1. Introduction. The classical Cramer-Lundberg model in insurance as- 
sumes that the risk process {R(t), t > 0) is given by R(t) = ct — Yi, 
where c> is the premium rate, (1^, i£N) is an i.i.d. sequence of nonneg- 
, ative random variables modelling individual claims, and (N(t), t > 0) is a 

homogeneous Poisson process of rate A > 0, independent of (Yj,, i £ N). Hence 
the cumulative claim process is modelled by the compound Poisson process 
J2i^i Yi- Let F denote the distribution function of Yi, and let fi = EY^. The 
central question for the model is the computation of the ruin probability in 
infinite time, given initial capital x > 0, defined by 

:= F(R(t) + x < for some t > 0). 

In case c < Xfi, this quantity is identically equal to 1. Hence, one always 
■ assumes the net profit condition c > A/U, and defines the parameter p := 
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Xfi/c < 1. Instead of studying the ruin probability, one can equivalently 
consider the survival probability 9{x) := 1 — which is more convenient. 
One of the few explicit results for the survival probability is the Pollaczek- 
Hinchin formula: 

oo 

(1.1) 9(x) = (l-p)^p n F?*(x), 

n=0 

where Fj(x) = (l/p) Jjf(l — F(t))dt is the integrated tail distribution. For- 
mula (1.1) is usually derived via renewal arguments. The resulting integro- 
differential equation for $ is solved using Laplace transforms. The expla- 
nation of the formula is given by considering the supremum of the dual 
process R(t) := —R(t). By the net profit condition R(t) drifts to — oo, hence 
the supremum is a.s. finite, and clearly, 6(x) = P(sup 0<i<oo R(t) < x). It is 
easy to see that sup 0<i<oo R(t) is a sum of geometrically many i.i.d. random 
variables. It is not, however, quite as easy to determine the distribution of 
these variables. Usually fluctuation theory is used. We refer the reader to [1] 
and [9] for details. 

In this paper we are interested in generalizations of the Cramer-Lundberg 
model, which lead to the same type of the Pollaczek-Hinchin formula for 
the survival probability, and which admit an explanation of the formula 
by decomposition of the supremum of the dual process in the random sum 
of ladder heights. One possible generalization of the model is to allow for 
additional uncertainties in the cumulative claims and/or in the premium 
income. These uncertainties may be the result of fluctuations in the claim 
arrival intensity, the number of insurees, inflation or surplus investment (see 
[9], page 568). Dufresne and Gerber [3] considered the risk process (R(t), t > 
0) perturbed by a multiple of standard Brownian motion (W(t), t>0), and 
defined X(t) := R(t) +sW(t), ? > 0. Using renewal arguments, they derived 
the formula 

oo 

(1.2) 6(x) = (l-p)J2 p n {G ( - n+1 > * Ff*)(x). 

n=0 

The parameter p and the distribution function Fj are the same as in the 
unperturbed model, while G is an exponential distribution function with pa- 
rameter 2c/? 2 . They also gave the following interpretation of the formula (1.2): 
Let 0i, 02) • ■ • be the moments when a new supremum of the dual process 
X(t) := —X(t) is reached by a jump of the claim process 

Eii )y i- Then the 

number of such moments has geometric distribution with parameter p, G is 
the distribution function of the supremum of X(t) just before a±, and Fj is 
the conditional distribution of the overshoot over the previous supremum, 
given 0i < oo. 
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Furrer [5] considered the process X(t) = R(t) + Z a (t), where R is the 
classical risk process, and Z a is an a-stable Levy process with no posi- 
tive jumps, 1 < a < 2. He used the explicit formula for the Laplace expo- 
nent of the infimum of X(t) due to Zolotarev [13] to obtain formula (1.2) 
for the survival probability of X(t). The distribution function G is ex- 
plicitly identified as the Mittag-Leffier distribution given by 1 — G{x) = 
YZ=v{-cx«- l Y/T{l + { a -l)n). 

Schmidli [11] gives a nice interpretation of G as the distribution of the 
supremum of the dual process X just before the first time the process X 
exceeds its previous supremum by a jump of the cumulative claim process. 
His setting is more general in the sense that the cumulative claim process is 
generated by a stationary, ergodic, marked point process. 

Another possible generalization of the classical risk process is to allow 
a different cumulative claim process. Dufresne, Gerber and Shiu [4] con- 
sidered the model in which the claim process was modelled by a Gamma 
process. Such a process has infinitely many jumps in finite intervals. Yang 
and Zhang [12] studied this model perturbed by a Brownian motion. Using 
the approach in [5], they derived a formula of the type (1.2) with Fj replaced 
by an exponential integral type distribution, and G is again the exponential 
distribution. 

The goal of this paper is twofold: (1) to extend the Pollaczek-Hinchin for- 
mula to the more general setting of the spectrally negative Levy processes 
while retaining the risk insurance theory interpretation; (2) to give a uni- 
fied and transparent approach to the problem by use of well-developed and 
powerful fluctuation theory for Levy processes. We will consider a general 
perturbed risk process X(t) = ct — C(t) + Z(t), where (C(t), t > 0) is a cu- 
mulative claim process, and (Z(t), t > 0) is a perturbation. Note that the 
cumulative claim process has to be increasing. Therefore, if one wants to 
stay in the realm of processes with stationary independent increments, the 
only choice for modelling (C(i), t > 0) is subordinators. Hence, we assume 
that (C(t), t > 0) is a subordinator (without drift) having finite expectation 
satisfying the net profit condition c — EC(1) > 0. The perturbation is mod- 
elled by a Levy process (Z(t), t>0) with no positive jumps, having zero 
expectation. The assumption that the expectation is zero is inconsequential, 
since EZ(1) can always be moved to the premium rate. In the analysis of 
the risk process (X(t), t>0), we will rely heavily on fluctuation theory for 
general Levy processes, which is particularly explicit for processes with no 
positive jumps. For background on these results, we refer the reader to the 
book by Bertoin [2]. 
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Our first result is the formula for the survival probability for the process 
X which is proved in Section 3: 



We essentially follow the approach from [5], and obtain explicitly the pa- 
rameter p and the distribution functions G and H appearing in the formula. 
It turns out that G can be identified as the distribution function of the 
absolute supremum of the process (—ct — Z(t), t > 0), while H is related 
to the subordinator C(t) only, and can be thought of as the integrated tail 
distribution of jumps. In Section 4 we give an interpretation of formula (1.3) 
by decomposing the supremum of the dual process X(t) := —X{t) into the 
random sum of modified ladder heights. In order to do this, we first show 
that the times when the new supremum of X(t) is reached by a jump of 
the subordinator are discrete. Let < o\ < 02 < ■ • ■ be those times, and 
let G be the distribution function of S(ai-), where S(t) := sup 0<s<t X (s) . 
We show that S(ai-) and the overshoot S(ai) — S(ai—) are conditionally 
independent given o\ < 00, and identify the conditional distribution of the 
overshoot with H. Using the strong Markov property at times o"j, we rederive 
formula (1.3) with G instead of G (and the same p). This clearly implies 
that G = G, yielding the required interpretation. Our results are more gen- 
eral and cover the results obtained in [5], [12] and [11] (in Levy case). 

Another interpretation of formula (1.3) is provided by looking at the lad- 
der height process of X. The ladder height process is obtained by time- 
changing S(t) by the inverse local time at zero of the reflected process 
S(t) — X(t). This process records only values where the new supremum 
is reached, and consequently, contains all the relevant information on the 
distribution of the supremum of X(t). In Section 5 the results of Section 4 
are reinterpreted and improved in terms of the ladder height process. 

We end this Introduction by noting that in a very recent paper Kliippel- 
berg, Kyprianou and Mailer [7] study ruin probabilities for general Levy 
insurance risk process (not necessarily spectrally negative) drifting to —00. 
They are mostly concerned with the asymptotic results for the first passage 
time and overshoot behavior at high levels. 

2. Setting and notation. Let (S1,JF, P) be a probability space on which 
all random variables will be defined. As explained in the Introduction, we 
model the cumulative claim process by a subordinator C = (C(t), t > 0) 
without a drift. Let v be the Levy measure of C; that is, v is a cr-finite 



(1.3) 



6(x) :=P inf X(t) > -x 




00 



= (1 - p) p n {G i - n+1) * * H n *)(x). 



n=0 
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measure on (0, oo) satisfying Jj ^(x A l)u(dx) < oo. The Laplace exponent 
of C is defined by 



so that 



Note that 



M/3):=/ (1-e-^Mds) 
J (0,00) 



E[exp{-pC(t)}] =exp{-t$ c (/3)}. 



/• /"OO 

EC(l) = $^(0+)= / xi/(ds)= / v(x,oo)dx, 
J (0,oo) Jo 

where the last equality follows by integration by parts. As explained in the 
Introduction, we assume throughout that EC(1) < 00. Let 

(2.1) daif"** ■»)*• 

Then if is an absolutely continuous distribution function with density ft,(x) = 
oo)/EC(l). We call the integrated tail distribution. The Laplace 
transform of H is given by 



/•oo /*oo 

(2.2) CH{(5):= e~ Px H{dx) = e -^ x h(x)dx 
Jo Jo 



1 <Ml) 

EC(1) /3 



Let AC(t) = C(i) - C(i-). It is well known that (AC(t) : t > 0) is a Poisson 
point process with characteristic measure z/ and state space (0, 00) U {d}. 
The cemetery state d is added to keep with the standard definition of a 
Poisson point process (cf. [8], page 435). It is assumed that the process is 
in state d whenever there is no jump and v({d}) = 0. Moreover, one has 

c(t) = j2 0<s < t AC(s). 

We model the risk process R = (R(t), t > 0) as R(t) = ct — C(t), where 
c > is the premium rate. Then R is a Levy process with no positive jumps 
(i.e., spectrally negative Levy process). The Laplace exponent ipR of R is 
defined by relation 

E[exp{/3i^)}]=exp{^ R (/3)}. 

Clearly, ipn((3) = c(3 — $>c(0)- It is important to note that R(t) stays positive 
in a neighborhood of t = 0, implying that ruin (with zero initial capital) does 
not occur immediately. This follows from the fact that lim^o C(t)/t = 
([2], page 84 or 192). From now on we assume that the net profit condition 
c > EC(1) holds, and let d:=c- EC(1). It follows that ER(1) = ^(0+) = 
d > 0, which implies that R drifts to +00. We also introduce the parameter 
p:=EC(l)/c 6 (0,1). 
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The perturbation Z = (Z(t), t > 0) of the risk process R will be modelled 
by a spectrally negative, mean zero, Levy process. Its Levy measure Hz is 
an infinite cr-finite measure on (— oo,0) satisfying the usual condition 

(2.3) / (x 2 A l)U z (dx) < oo, 

J(-oo,0) 

and the additional condition 

(2.4) / \x\Ii z {dx) < oo, 

J(-oo,-l) 

which ensures finite expectation of Z. The Laplace exponent of Z is given 
by 

(2.5) ^z(P) := ^/3 2 + / (eP* - 1 - f3x)Yl z (dx), 

2 J(-oo,0) 

where c > 0, and integrability of the integrand follows from condition (2.4). 
Further, EZ(1) =ip' z (0+) = (e.g., [10], page 163). Note that we allow Z 
to be identically zero (both 11^ = and c = 0). However, Z cannot be com- 
pound Poisson because such processes cannot have KZ(t) = 0. Let us point 
out that our setting includes the Brownian perturbation (c > 0, Hz = 0), 
and also the perturbation by a-stable spectrally negative Levy process for 
a G (1, 2) (c = 0, U z (dx) = (a/|x| Q+1 )l ( _ OOi0 ) dx). 

Finally, we define the general perturbed risk process X = (X(t), t>0) as 

X(t) := R(t) + Z(t) =ct- C(t) + Z(t), 

where C and Z are independent processes. The process X is a spectrally 
negative Levy process with finite positive expectation EX(1) = c — EC(1) = 
d > 0. Therefore, lim t ^ 00 X(t) = +oo a.s., that is, X drifts to infinity. Let 
^0(t) :=a{C(s),Z(s), < s < t), and let T = (F(t), t > 0) be the filtration 
obtained in the usual way by augmenting J^(t). Clearly, X{t) is J r (t)-measurable. 
The Laplace exponent ip of X , defined by the relation 

E[exp{/3X(t)}]=exp{tV(/5)}, 
is, due to independence of C and Z, given by 

il>(J3)=cP-$c(fi)+il>z(P), P>0. 

Since ip is strictly convex and t//(0+) = EX(1) > 0, ip is strictly increasing 
on [0, oo), and therefore has a strictly increasing inverse $ : [0, oo) — > [0, oo). 
Since ^(0) = 0, it follows that $(0) = 0. 

In the sequel, we will be interested in the function 8 : [0, oo) — > [0, 1] defined 

by 



(2.6) 



d{x) := F(X(t) > -x, for all t > 0). 
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This function is the survival probability of the general perturbed risk process 
X starting with the initial capital x > 0. The initial behavior of X determines 
9 at zero. If there is no perturbation, that is, if X = R, then, as said before, 
X remains positive (a.s.) for an initial period of time, and hence 6(0) > 0. On 
the other hand, if Z ^ 0, then X is of unbounded variation, hence the point 
is regular for (— oo,0) ([2], page 192). Thus X hits the interval (— oo,0) 
immediately, implying 6(0) = 0. 

3. Laplace transform approach. In this section we derive the Pollaczek- 
Hinchin formula for the survival probability using the explicit form of the 
Laplace transform of the absolute infimum of X. Let I(oo) := info< s <oo X(s) 
and I(t) := info< s <t X(s). The fluctuation theory for Levy processes pro- 
vides the following formula for the Laplace transform of the infimum evalu- 
ated at an independent exponential time r(q) with parameter q > (see [2], 
page 192): 

E ww^)))H t g;g, r $>o. 

F 

Letting q { 0, and using I(r(q)) — > /(oo), it follows that 

(3.1) E[exp{/?I(oo)}] = V>'(0+)^y = d-^ , (3 > 0. 

Let us introduce for a moment the following notation: Y(t) = ct + Z(t) and 
ipY (P) = c/3 + ipz (/3) • By the same argument as above it follows that 



(3.2) E 



Let G denote the distribution function of — info<i <0 o Y(t) = sup 0<<<00 (— ct- 
Z(t)). Then the last formula says that 

(3.3) £G(f3):= e -f 3 *G(dx) = c-P— f3 > 0. 

./o W(p) 

Recall formulae (2.1) and (2.2) from Section 2: 

H(x):=(l/EC(l)) [\(y,oc)dy and CH(f3) = $ c (/?)/(EC(l)/3). 



o 



Also recall that p = EC(l)/c, hence d/c = (c - EC(l))/c = 1 - p. Now we 
compute d(3/ip((3) in terms of p, CG and CH . This idea comes from [5]: 

d4^= d 



1 



d 



c/CG(p) -EC(l)CH(p) 
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_ d £G(p) 

~ c 1 - P CG{(3)CH{(3) 

oo 

= (l-p)JM(J})Y,(pCG(J3)CH(J3)) n - 

n=0 

By inverting the Laplace transform, we obtain the following theorem. 

Theorem 3.1. The survival probability of the general perturbed risk pro- 
cess X is given by 

9{x) = P(/(oo) > -x) 

(3.4) 



(1- p)Y / p n (G (n+1) * *H n *){x), x>0. 



n=0 

We point out that H depends only on the subordinator C, while G de- 
pends on the premium rate c and the perturbation Z. Brownian perturba- 
tions were considered in [5] and a-stable ones in [3] and [12]. In both cases 
the distribution G is given explicitly. If there is no perturbation, Z = 0, then 
CG{(3) = 1, and consequently, the distribution function G can be omitted 
from formula (3.4). 

4. Decomposition of the supremum of X. Let X(t) := —X(t) = —ct + 
C{t) — Z{t) denote the dual process of X. Let 

S(t) := sup X(s) and S(oo) := sup X(s). 

0<s<t 0<s<oo 

Since X drifts to — oo, S(oo) < oo a.s. Introduce the following notation: 
I(t) := inf < s < t X(s) and S(t) := sup < s < t X(s). Clearly, -I(t) = S(t). By a 

time reversal argument, —I(t) = S(t) — X(t), and hence 
(4.1) S(t)-X(t)±S(t). 

In this section we give a decomposition of X at certain stopping times which, 
following Schmidli [11], we call modified ladder epochs. 

Let V(J-) be the predictable cr-algebra on K + x with respect to the 
filtration T introduced in Section 2. Let Bq be the Borel a-algebra on 
(0, oo) U {d}. If H : R + x Q x ((0, oo) U {d}) — > R + is a nonnegative process 
measurable with respect to V(J-) ® Bq, then the following compensation 
formula is valid (e.g., [8], page 439, or [2], page 9): 



E 

(4.2) 



( H(t,co,AC(t,cu))\ 

\0<t<oo / 

= E ( f°° dt f v{de)TL{t, u, e)] . 

\Jo J(0,oo) J 
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The first use of this formula will be to compute the expected number 
of times the new supremum of X is attained by a jump of a subordinator 
C over the previous supremum. Note that this is the case if and only if 
AC(t)>S(t-)-X(t-). 

Theorem 4.1. The following formula is valid: 
(4-3, e( Y. 1 i ^ " EC(1> 



\0<t<oo 



L {AC(t)>S(t-)-X(t-)} J c _ EC(1) ' 



Proof. Take H(t,LO,e) := ^rgn (t- w) oo)^ * n ^ e com P ensa tion 

formula. The left-hand side in (4.2) is then precisely the left-hand side 
in (4.3). For the right-hand side in the compensation formula, compute 

E (r dt /o,oo) Kde)i ^-)-^-)-) (e) 



E (^°° dtv(S(t-)-X(t-),oo) 
E[v(S{t) - X(t),oo)]dt 



(4.4) 



diE[i/(S(t),oo)], 

where the third line follows by continuity in probability of X , and the fourth 
line by (4.1). Clearly, the last expression is equal to the monotone limit 



(4.5) 



roc 

lim / e- qt dtE{v(S(t),oo)] 



1 f°° 

lim - / qe~ qt dtE[u(S(t),oo)}. 
q^o q Jo 

Let r(g) be an exponential time with parameter q independent of C and Z, 
and let F denote the distribution function of S(r(q)). Then F is exponential 
with parameter &(q). It follows that 

qe- qt dtE[u(S(t), oo)} 

qe- qt dtu(S(t),oo) 



= E 

(4.6) = E[i/(S(t(?)),oo)] 

u(x, oo)F(dx) 

' (l-e~^ x )u(dx), 
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where the last equation follows by integration by parts. Further, 

l_ e -*(?)« $r q \ i x 
hm = lim x = — — -x = — . 

q^o q q->o q ^ (0+) d 

By the monotone convergence theorem, 

POO 

lim / e - qt dtE[u(S(t),oo)] 



(4.7) 



/■°o ^ _ e -*(g)s 

lim / 



xv{dx) 
EC(1) 



c-EC(l) 
This proves formula (4.3). □ 

Remark 4.2. We would like to emphasize a very important and some- 
what subtle point which is a consequence of Theorem 4.1. Namely, the epochs 
when a new supremum of X is reached by a jump of C are discrete, and, 
in particular, neither time nor any other time is an accumulation point of 
those epochs. More precisely, let us define 

(4.8) en = a := inf {t > : AC{t) > S(t-) - X(t-)}, 
and inductively, 

(4.9) a n+l := inf {t > a n : AC(t) > S(t-) - X(t-)} 

on {<7 n < oo}. Theorem 4.1 implies that a > a.s. and a n < <r n +i a.s. on 
{<j n < oo}. As a consequence, we can order the epochs when a new supre- 
mum is reached by a jump of a subordinator as < o\ < 02 < • • • a.s. The 
decomposition of X and the ensuing derivations will depend on this result 
in a crucial way. 

For y > 0, let r y := inf{i > 0:X(t) > y} be the entrance time of X in 
(y,oo), and, similarly, r y := inf{t > 0:X(t) > y}. Note that S(t—) < y if 
and only if t<f y . We need the expected occupation time formula for the 
reflected process 5 — X before a Ar y . 

Proposition 4.3. For x > and y > 0, the following formula is valid: 

(4.10) e( ATs \s ( t)-x mx) dt = F ^ = ^ ? y = °°) J 
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Proof. We first compute the expected occupation time of S — X below 



x: 



E 



1 (S(t)-X(t)<x) dt 



F(S(t)-X(t)<x)dt 



(4.11) 



o 



¥(S(t)<x)dt 



Since (r x , x > 0) is a subordinator with the Laplace exponent it follows 
that Et x = (Enjx = &(0+)x = x/d. 

Now we compute the expected occupation time of S — X below x after 



time a V r y : 



roo 

E / 1 



(%-X(t)< :C ) 1 (*>^) 1 (5(t)>y) dt 



E 



1 



a\/r v 



(S(t)-X(t)<x) 



dt, a V T y < oo 



(4.12) 



(er V f y < co)E 



dt\a V f„ < oo 



aVr v 



= P(a Vr y < oo)E ^ l^^j dt 

X 

= P(cr < oo, < °°)^- 

To justify the passage from the third to the fourth line, note that a V f y is a 
stopping time at which X reaches a new maximum, and hence by the strong 
Markov property, the reflected process S — X starts afresh from 0. Similarly, 



(4.13) 



E / Q 1 (S(t)-X(t)<x) 1 (S(t)>y) dt 



V 

Ty <00)~. 



Subtracting (4.12) from (4.13), it follows that 



(5(i)-X(t)<x) 1 (t<^) 1 (5(i)>j/) 



roo 

(4.14) E / 1 
Jo 

One can prove similarly that 
(4.15) 



.r 



dt = P(<7 = OO, T v < OX)) — . 

y d 



roo 

E / t, 
Jo 



,t (t<a) dt = ¥(a = oo)- 



L (S(t)-X(t)<x)Ht<*)'~- -v- 
Finally, (4.10) follows by subtracting (4.14) from (4.15). □ 



Note that the proposition says that the expected occupation time mea- 
sure for S — X before a A T y is proportional to the Lebesgue measure on 
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[0,oo). Hence, formula (4.10) is by definition of the expected occupation 
time measure equivalent to 

(4.16) E r ATy f( S(t) - X(t)) dt = F((J = 00 fy = °° ) r f{u)du, 
Jo a Jo 

where / is a nonnegative Borel function on [0, 00). 

Let J := (AC (a) — (S(a—) — X(a—)))tr a<QO \ be the overshoot at time a. 
In the next proposition we compute the preliminary version of the joint 
distribution of the vector (S(a—),J,S(a—) — X(a— )) on {a < 00}. 

Proposition 4.4. For x,y,z> 0, 

P(S(cr-) <y,J>x, S(a-) - X(a-) >z,a<oo) 

( 4 - 17 ) p , ~ , 

F(a = 00, T y = 00) 



d 



00 

u(u, 00) du. 

x+z 



Proof. We use the compensation formula with 
H(t,u,e) 

Then 

0<i<oo 

= P(5(<t-) < y, S(a-) - X(a-) > z, J > x, a < 00). 
On the other hand, 

E[ [°° dt [ v(de)H(t,u,E)) 

\Jo 7(0,00) / 



' / dtl (S(t-)<y) 1 (S(t-yx(t^)>z) 



(0iOo) 1 (, + 5(^)-X(^),oo)( £ )^ & : 

E Uo ATV dt Hsv-xit^fi* + 8® ~ X(t), OO)) 

\z,oo)( u ) v ( x + u i 00 ) du 
v(u, 00) du, 



d 

(a = 00, T y = 00) 



d 



X + Z 
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where the fourth line follows from (4.16) with f(u) = lr z 00 )(u)^(x + u, oo). 
□ 

From formula (4.17) we can easily derive several useful corollaries. 

Corollary 4.5. The following formulae are valid: 

(4.18) P(cr<oo) = p, 

l r 00 

(4.19) P( J > x\a < oo) = / i/(u, oo) du = 1 - if (x). 

EG (1) Jx 

Proof. Let x — » 0, w — > oo and z — > in (4.17). It follows that 
P(.<oo) = IP((J J 00) EC(l). 

Solving for P(a < oo), we get (4.15). To obtain (4.19), let y — > oo and z — > 
in (4.17). It follows that 



, P (T = oo r , 

'( J > x, a < oo) = / u(u,oo)du. 

d Jx 



By conditioning, 



1 — p roc 

P(J > x\a < oo) = — / vyu, oo) du 

pd J x 



1 



EC(1) 



f(u, oo) du. 



□ 



In the next corollary, we interpret S(o~—) as the absolute supremum 5(oo) 
in case a = oo. 

Corollary 4.6. The event {a < oo} and the random variable S(o~—) 
are independent. As a consequence, the conditional distribution of S(cr—) 
given a < oo is equal to the unconditional distribution of S(o~—). 

PROOF. Let x — > and z — ► in (4.17). It follows that 
P(5(<r-) <y,cr<oo) 

EC(1) 



(4.20) =p(cr = oo,f y = oo) 



:P(<r = oo,S(oo) <y) 



d 

EC(1) 



d 

Clearly, 

P(5(<r-) <u,ct = oo) =P(5(oo) <y,<r = oo). 
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Adding up, 

P( S(a-) <y)= + l) P( S(oo) < y, a = oo) 

= ^F(S(oo)<y,a = oo). 
a 

Therefore, 

P(5(<t-) < y)P(<7 < oo) = ^P(S(oo) < y,a = oo)^^ 

= P(S(o— ) <y,a<oo) 

by (4.20). □ 

It follows that 

P(cr = oo, T y = oo) = P(er = oo, 5(oo) < y) 
= F(a = oo,S(a-)<y) 

(4-21) 

= P(a = oo)P( S{a-) < y\a = oo) 
= P(<t = oo)P(S(<t-) <y). 

Let G denote the distribution function of S(a—). Proposition 4.4 can be 
now improved as: 

Theorem 4.7. The distribution of the vector (S(a—), J, S(a—) — X(a—)) 
on the set {a < oo} is given by 

P( S(a-) <y,J>x, S(a-) - X(a-) > z, a < oo) 

(4.22) 

= P( 5(<7-) < y) (^^y !/(u, oo) duj P(a < oo) . 

Moreover, S(a-) and J are conditionally independent given a < oo, and 
(4.23) P( S(a-) < y, J > x\a < oo) = G{y){l - H{x)). 

Remark 4.8. Formula (4.22) considerably extends the severity of ruin 
formula (see, e.g., [9], page 168). 

It is now possible to write the absolute maximum of A" as a random sum 
of modified ladder heights. Recall that o~\=o and 

a n+1 = mf{t > o n : AC(t) > S{t-) - X(t-)} 
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Fig. 1. A trajectory of the dual process X(t) = —ct + C(t) + <;W(i), where C(t) is a 
compound Poisson process, W(t) is a standard Brownian motion and ? > 0. 

on {a n < oo}. Let Lq := 5(oi— ), J\ := S{o\) — S{a\—) and L\ := S{o2— ) — 
S(ai) on {a i < oo}, and so on. We call Lq, Ji, Li, . . . the modified ladder 
heights. Let also iV := max{n:cr n < oo}. By strong Markov property of X, 
N has geometric distribution with parameter ¥(a± = oo) = 1 — p. Clearly, 

(4.24) S(oo) = L + Ji + Lx + • • • + J N + L N . 
See Figure 1. 

Note that P(L < x,N = 0) = P(5(a-) < x,cj = oo) = G(x)(l - p). For 
every n € N, by the strong Markov property at a n , and by (4.23), we have 

P(L + Ji+L l + --- + J n + L n <x,N = n) 
= {l- p)p n {& n+l) * *H n *)(x). 

This leads to the Pollaczek-Hinchin formula for the distribution function 
of S(oo). 

Theorem 4.9. For x > 0, 

oo 

(4.25) P( S(oo) <x) = (l-p)J2 P n {& n+l) * * H n *)(x). 

n=Q 
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Equating this formula for 9{x) = P(5(oo) < x) and (3.4), one obtains 

oo oo 

(1 - p) p n (G ( - n+1 > * H n *)(x) = (l-p)J2 p n (& n+1) * * H n *)(x). 

n=0 n=0 

By computing Laplace transforms of both sides, we get that 
(496) (1-pKG(/3) _ (I-P)CG(P) 

from which it immediately follows that G = G. Thus we have proved the 
following 

Corollary 4.10. The random variables 

sup {-ct-Z{t)) and sup {-ct + C(t) - Z(t)) 

0<t<oo 0<t<a 

have equal distributions. 

5. Decomposition of the ladder height process. In the previous section 
we looked at the process X at^ the modified ladder epochs o\ < o-i < ■ ■ ■ , and 
essentially decomposed the X at these epochs. In this section we consider 
the ladder height process if of A" obtained by time-changing the supremum 
process S by the inverse local time at zero of the reflected process S — X. 
The excursion representation of the process S — X will be combined with 
fluctuation identities and results from Section 4 to give a detailed description 
of the ladder height process. _ _ 

Let us first briefly recall that S — X is a strong Markov process, and 
hence it admits a local time process at 0, L = (L(t), t > 0). The process 
L is continuous and nondecreasing, and increases only when S — X is at 
zero, or in other words, when X reaches a new supremum. More precisely, 
the support of the measure L{dt,uj) is the zero set of S(t,u>) — X{t,oj). If 
L -1 (t) := inf{s > : L(s) > t} is the inverse of L and one defines 

S(L~ 1 (t)), L" 1 (t)<oo, 
+oo, otherwise, 

it is well known that the process {{L~ l (t) , H (t)) , t < L(oo)) is a two-dimensional 
subordinator killed at rate q := l/EL(oo) ([2], page 156). In particular, 
H = (H(t), < t < L(oo)) is a subordinator killed at rate q = l/EL(oo). 
Clearly, 

(5.2) H(L(oo)-)= sup H(t)= sup S(t) = S(oo), 

0<t<L(oo) 0<t<co 

and hence, the distribution function of S(oo) is equal to the distribution 
function of H(L(oo)—). Fluctuation identities give a formula for the Laplace 
exponent of H . 



(5.1) H(t) := | 
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Lemma 5.1. The Laplace exponent k of H = (H(t), < t < L(oo)), with 
L suitably normalized, is given by the following formula: 

(5.3) k(P) = ^ =d + EC(l) / (l-e- f3x )H{dx) + ^f^, 

P 7(0,oo) p 

where H is the finite measure defined in (2.1). 

PROOF. The bivariate Laplace exponent /?(a,/3) of {{L~ l (t) , H (t)) , t < 
L(oo)) is defined by 

exp{-£(a,/3)} = E[exp-{a£ _1 (l) + /3H(1)}], a,/3>0. 
The explicit formula for k comes from fluctuation theory: 

a - </>(/?) 



(5.4) K(a,f3)=k 



where fc is a constant depending on the normalization of the local time. We 
take k = l. By letting a = in (5.4), we obtain the Laplace exponent of H: 



(5.5) 



£(/?) = k(0,/?) = ^M 



/3 

Further, integrating by parts, we get 



+ ^^, /3>0. 



c+ / (1 — e ^ x )i^(x, oo) dx — / z^(a;,oo)(ia; 



roo 

= c + EC(l) / (l-e- /3:c )i/(dx)-EC(l). 
Jo 

Together with (5.5), this gives (5.3). □ 

Remark 5.2. Note that the same argument shows that the ladder height 
process of (—Z(t), t>0) has the Laplace exponent equal to ipz(P)/P- 

Let L be the local time of S — X normalized by the choice k = 1 in for- 
mula (5.4). The excursion process (e s :s > 0) of the reflected process S — X 
can be viewed as the superposition of three independent Poisson point pro- 
cesses: finite duration excursions that end with a jump of the subordinator, 
finite duration excursions that do not end with a jump of the subordinator, 
and excursions of infinite duration. Note that one can include the jump, if 
any, that concludes an excursion as part of that excursion and retain a Pois- 
son point process. Also one needs to "carry along" the information about 
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which jumps come from C and which from Z, but that is a question of 
choosing a suitable filtration. 

The excursion process and with it the ladder height process is killed at 
the time of the arrival of the first excursion of infinite duration. From the 
excursion picture we know that L(oo) is an exponential random variable 
with parameter equal to the killing rate of the ladder height process. By 
convention, H(t) = oo after killing. It is easily shown that the killing rate 
can be obtained from the Laplace exponent as ■0(0+). Since EZ(t) = 0, we 
know that ipz(P) /P ~ ¥ as (3 — > 0. The measure H is finite, so by dominated 
convergence, 

/ (1 - e-P x )H(dx) -► 

i(0,oo) 

as — ► 0. Using the explicit formula (5.3) one finds that the killing rate 
equals d = c — EC(1), or in other words, L(oo) ~ exp(c — EC(1)). 

Recall that a is the first modified ladder epoch. On the local time scale 
L(a) corresponds to the time of the first arrival of an excursion that ends 
with a jump of the subordinator C unless the excursion process is killed first. 
The probability ¥(a < oo) is therefore equal to the probability that the first 
excursion that ends with a jump of the subordinator arrives before the first 
excursion of infinite duration. We are thus computing the probability that 
one of the two independent Poisson processes will "claim" the first arrival. 
It is well known that we can compute 

p = P(<7 < oo) =¥{L(a) <L(oo)), 

where L(a) and L(oo) are two independent exponential random variables. 
Since we know the rate of L(co), a simple calculation shows that the rate of 
arrival of excursions that end with a jump of C equals EC(1). 

Let us turn to the ladder height process H. The jumps of H are a Poisson 
point process. These jumps, however, can be seen as a mapping of the "big- 
ger" process e of excursions. By the mapping theorem (see [6], page 17) the 
jumps of H are a superposition of two independent Poisson processes killed 
at an independent exponential time of rate d = c — EC(1). The mapping the- 
orem is applied in such a way that the image of e s is d if the excursion ends 
without a jump. The jumps coming from excursions that end with a jump 
of C contribute a pure jump process to H . The other excursions contribute 
jumps and possibly a deterministic drift. But the jumps coming from Z are 
an independent Poisson process. Hence the process H is a sum of two in- 
dependent subordinators rj and £ killed at an independent exponential time 
r := L(oo). 

The subordinator ij corresponds to increases of S due to jumps of C. As 
the set of times when S increases by a jump of C is discrete, the process 
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rj is compound Poisson with arrival rate equal to EC(1). Jumps of r\ have 
the same distribution as the overshoot J which is given by the formula 
(4.19): P(J € dx\a < oo) = u(x, oo) dx/KC(l). This means that 77 contributes 
exactly 



to the Laplace exponent of the ladder height process H given in (5.3). 

The subordinator £ arising from increases of H not due to jumps of C is 
independent of 77. This leaves us with the conclusion, given that the killing 
rate is c — EC(1), that the Laplace exponent of Q is ipzift)/ P- This way the 
groundwork in Section 4 has been translated into a decomposition of the 
ladder height process: 

Theorem 5.3. Let X(t) = -ct + C{t) - Z(t) and let H be the ladder 
height process of X. The following assertions are valid: 

(i) H is killed at rate d = c — EC(1). 

(ii) The ladder height process is the sum of two independent subordinators 
(H(t) = r](t) +C(0 : < t < L(oo)). The subordinator rj corresponds to jumps 
of S due to the claim process C . It is compound Poisson with Levy measure 
KC(1)H. The subordinator £ corresponds to increases of S not arising from 
jumps of C. Its Laplace exponent is ipz(P)/fl- 

The decomposition of H gives some insights into the structure of the 
process H. From the form of the Laplace exponent of we find that the 
distribution of £ does not depend in any way on the distribution of C, which 
is a remarkable conclusion given that C is a subordinator with a dense set 
of times of jumps. Many other conclusions from Section 4 can be recast in 
terms of the ladder height process. 

As an example, consider Corollary 4.6. Recall that r = L(oo). Define 7 := 
L(a) if a < 00 and 00 else. The event {a < 00} is equal to {7 < r}. It is 
easily seen that 



is G. The decomposition of H gives further information about G. Since the 
random variables 7 and r are independent of (, and 7 A r ~ exp(c), the 
distribution of S(a—) is that of a subordinator taken at an independent 
exponential time and hence infinitely divisible with Levy measure A given 




(0,00) 
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by 

(5.6) A(dx)= [ - — P(CtGds)dt. 

J (0,06) t 

See [2], page 162. However, a direct computation using the Laplace exponent 
of C yields 

(5.7) E(exp(-/?C 7 Ar))= / E(exp(-/3Ct))ce- rf ^ 

J(0,oo) 

= c[ e^p(-tip z {f3)/f3)e~ ct dt 

J(0,oo) 

Comparing this formula to (3.2) gives an independent proof of Corollary 4.10. 
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